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The adsorption of a bottle-brush polymer end-grafted with one chain end of its backbone to a flat 
substrate surface is studied by Monte Carlo simulation of a coarse-grained model, that previously has 
been characterized in the bulk, assuming a dilute solution under good solvent conditions. Applying 
the bond fluctuation model on the simple cubic lattice, we vary the backbone chain length Nb from 
Nb — 67 to Nb = 259 effective monomeric units, the side chain length N from N = 6 to N = 48, and 
the grafting density a = 1, i.e., parameters that correspond well to the experimentally accessible 
range. When the adsorption energy strength e is varied, we find that the adsorption transition (which 
becomes well-defined in the limit Nb — > oo, for arbitrary finite N) roughly occurs at the same value e c 
as for ordinary linear chains (JV = 0), at least within our statistical errors. Mean square end-to-end 
distances and gyration radii of the side chains are obtained, as well as the monomer density profile 
in the direction perpendicular to the adsorbing surface. We show that for longer side chains the 
adsorption of bottle-brushes is a two step process, the decrease of the perpendicular linear dimension 
of side chains with adsorption energy strength can even be non-monotonic. Also the behavior of the 
static structure factor S(q) is analyzed, evidence for a quasi-two-dimensional scaling is presented, 
and consequences for the interpretation of experiments are discussed. 



I. INTRODUCTION 

Macromolecules with a comb-like chemical architec- 
ture, where flexible side chains are densely grafted to a 
(also flexible) polymer acting as "backbone" have found 
abiding current interest, see, e.g., the recent reviews 0- 
4]. These so-called "bottle-brush polymers" exhibit an 
interesting competition due to the steric repulsion be- 
tween the side chains and the configurational entropy of 
the backbone: varying the grafting density of the side 
chains and their length, the effective stiffness of these 
cylindrical brushes can be controlled over a wide range. 
Since the overall conformation of the bottle-brush poly- 
mer is sensitive to external stimuli such as solvent quality, 
pH value and ionic strength of the solution, or tempera- 
ture, electromagnetic fields, various possible applications 
of these bottle-brush polymers have been discussed (ac- 
tuators, sensors, building blocks of new nanostructures or 
templates for producing metallic nano- wires. 0-0). A 
fascinating aspect is also the importance of biomolecules 
with bottle-brush architecture, such as proteoglycans 
(polyelectrolytes consisting of a protein backbone and 
carbohydrate side chains, performing biological functions 
from cell signaling and cell surface protection to joint lu- 
brication Hi El). Of course, under many circumstances 
these functions of bottle-brushes occur when they are 
attached to a substrate surface (e.g. a cell membrane 
in a biological context, or an inorganic flat solid surface 
for nano-technological applications, including also special 
surface coatings [13]). Bottle-brush molecules attached 
to surfaces are also of particular interest, since additional 
experimental tools become available to study their struc- 
ture (e.g., one can directly visualize their large-scale con- 
formation by scanning force microscopy [H, EH Ell; or 



one can use atomic force microscopes to measure force 
vs. extension curves [l6j]). Thus, while rich information 
on the structure of bottle-brush polymers at surfaces al- 
ready is available @, it is not so clear which of these 
features are intrinsic properties of these complex macro- 
molecules, and which features are only induced by their 
adsorption to the substrate. We recall that basic aspects 
of bottle-brush polymers in dilute solution away from 
any constraining surfaces are still incompletely under- 
stood (e.g. the relation between the persistence length of 
these macromolecules and their microscopic characteris- 
tics such as side-chain and backbone chain lengths N and 
N b has been controversial, see e.g. [I3-E1)- There are 
already some indications from theoretical work [2(J [H[ 
that strongly adsorbed bottle-brushes exhibit properties 
different from the bulk. 

Thus the present paper intends to contribute to a clar- 
ification of this problem, by presenting a comprehen- 
sive computer simulation study of (large) bottle-brush 
polymers near surfaces, over the full range extending 
from a repulsive surface (applying only the constraint 
that one backbone chain end is grafted to this surface) 
through the region of the adsorption transition, where 
the overall conformation of the macromolecule changes 
from a three-dimensional "mushroom" to a quasi-two- 
dimensional "pancake" , to the strongly adsorbed almost 
two-dimensional case. Using side chain lengths up to 
N = 4:8 we ensure that statistical concepts of theoretical 
polymer physics HI start to become applicable HI , and 
we also note that our range of side chain lengths can be 
nicely mapped to experiments (l9j . While the adsorption 
transition of flexible linear chains has been studied exten- 
sively by various theoretical methods [2^ - l32| , the present 
work is the first one to present large-scale simulations on 
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the adsorption of bottle-brush polymers. Clearly, the in- 
terplay of the conformational changes due to adsorption 
with the chain stiffness (tunable via the side chain length 
variation) is of particular interest (33j . 

The plan of this paper is as follows: Section II describes 
the model and the simulation techniques while Section 
III presents the numerical results. Section IV gives a 
summary and an outlook on both pertinent experimental 
work and on open questions. 



II. MODEL AND SIMULATION TECHNIQUE 

The bond fluctuation model on the simple cubic lattice 
[33 - l37| has been extended in previous work to simulate 
bottle-brush polymers in dilute solution under good sol- 
vent conditions [19l. |38||. In the present paper, we extend 
these studies to bottle-brush polymers tethered with one 
chain end of their backbone to a flat impenetrable sur- 
face, which is chosen to be the xy-plane (z = 0) of the lat- 
tice. Besides the excluded volume interactions between 
the monomers of the side chains and/or the backbone 
chain, and the infinitely repulsive interaction between 
monomers and the grafting surfaces (which prevents that 
monomers can occur with z-coordinates z < in the 
system) , we also consider the effect of an additional at- 
tractive short range interaction e between the grafting 
surface and the monomers of the bottle-brush polymer. 
Recall that in the bond fluctuation model each effective 
monomer blocks all 8 sites at the corners of an elemen- 
tary cube of the lattice from further occupation. Thus, 
an energy e is won if 4 sites of such a cube are in plane 
z = 0, while the remaining 4 sites of the cube then have 
to be in the adjacent plane z = 1 (note that we use the 
lattice spacing as the unit of length). 

The bond vectors connecting two adjacent monomers 
along a chain are chosen from the set {[2,0,0], [2,1,0], 
[2,1,1], [2,2,1], [3,0,0], and [3,1,0]}, including also all pos- 
sible rotations, reflections and reversions of these bond 
vectors. We use the same set of bond vectors, irrespec- 
tive of whether we consider side chains of the bottle-brush 
or the backbone chain. Note that this model is one of 
the standard models used for Monte Carlo simulation 
of macromolecular systems, many physical properties of 
polymers are rather well accounted for, and very efficient 
simulation algorithms can be formulated for this model 

The chemical architecture of a bottle-brush polymer is 
arranged such that at a backbone monomer that acts as 
a grafting site one side chain is anchored, and the chem- 
ical distance (i.e., difference in the consecutive labels of 
the backbone monomers) is the inverse of the grafting 
density, I /a. Thus, only integer numbers I /a are pos- 
sible, and side-chains are anchored at regular chemical 
distances. Note that both next to the first grafting site 
and next to the last grafting site of the backbone we add 
one extra backbone monomer. Thus if we have n c side 
chains, the total number of monomers in the backbone 



is Ni, — \(n c — l)/cr + 1] + 2, and the total number of 
monomers in the bottle-brush then is JVtot = -W& + n c N, 
when N denotes the number of monomers per side chain. 
Now, we also require that one of the two extra monomers 
at a chain end of the backbone is attached to the adsorb- 
ing surface at z = 0, and hence our simulation deals with 
a "bottle-brush mushroom" (recall that an isolated poly- 
mer chain |24| , end-grafted to a planar surface tradition- 
ally is called a "polymer mushroom" , and we extend this 
nomenclature to bottle- brushes). In the present work, 
only a = 1 is considered. 

The initial configuration of the bottle-brush polymer 
is constructed by assuming that the structures of the 
backbone chain and of the side chains are both rod- 
like. Namely, the backbone is oriented in the direction 
along the + z-direction, setting all bond vectors between 
monomers equal to r = 3z, where z is a unit vector along 
the z-axis. The bond vectors between monomers on each 
side chain are chosen randomly from one of the allowed 
bond vectors in the x-y plane. This initial configuration, 
of course, is very far from equilibrium, and needs to be 
carefully relaxed towards equilibrium in the first part of 
the Monte Carlo run. 

In order to speed up the dynamics of the model, we 
do not move the monomers from their previous posi- 
tions to only one of the six nearest neighbor sites of the 
monomer, as is traditionally done, but allow attempts to 
move a monomer to one of the 26 sites surrounding the 
current position of the monomer. This "L26 move" al- 
lows bonds to cross each other and has been shown to 
yield a significant speedup of the algorithm in compar- 
ison with the standard "L6" move 1391] . In addition to 
the local moves, the pivot algorithm |40| is employed in 
the following ways: (a) A monomer on the backbone is 
chosen randomly and the part of the bottle-brush poly- 
mer containing the free end of the backbone is exposed to 
a trial move, by randomly choosing one operation from 
the 48 symmetry operations of the lattice (no change, 
rotations by 90° and 180° angles, reflections, and inver- 
sions). Of course, the trial configuration is accepted as 
a new configuration only, if it does not violate any ex- 
cluded volume constraints and if it passes the Metropolis 
acceptance test (where the change in adsorption energy 
enters). The same is true for the attempted "L26" moves, 
of course, (b) A monomer is chosen randomly from the 
monomers on all side chains, and the trial move is ef- 
fected by transforming the part of the side chain from 
the selected monomer to the free end by one of the same 
48 symmetry operations as mentioned above. 

A Monte Carlo step (MCS) in our simulation then con- 
sists of a sequence of N to t "L26" moves (each monomer 
on average is attempted to be moved once), k p b pivot 
moves of choosing a backbone monomer randomly, and 
kp S pivot moves of choosing a side chain monomer ran- 
domly. We set k ps — n c /A and adjusted k p b such that 
the acceptance ratio was about 30 % to 40% for small 
values of e. However, when a bottle-brush polymer be- 
comes strongly adsorbed to the surface successful pivot 
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FIG. 1. Snapshots of typical configurations of bottle-brush polymers for backbone chain length Nb = 131 and side chain length 
N = 12, and four choices of adsorption energy e, namely e = (a), 1.0 (b), 1.5 (c) and 2.0 (d). 
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FIG. 2. (a) Fraction of monomer-surface contracts N s /Ntot and (b) ratio R^ z /R 2 g n of the mean square gyration radii of 
the whole bottle-brush polymers perpendicular and parallel to the surface plotted against the adsorption energy e. Data for 
backbone lengths Nb = 67, 131 and 259 are included, as well as side chain length N = 0, 6, 12 and 24. 



moves are more difficult. To avoid the problem that in 
the regime of the strong adsorption the bottle-brush con- 
figurations are frozen, 16 to 512 different equilibrium con- 
figurations are generated as initial configurations for the 
start of Monte Carlo simulations in equilibrium. During 
the equilibration process, the convergence of the time se- 
ries for the energy and the gyration radius components 
of the bottle-brush polymer are monitored. For the mea- 
surements at each parameter set in total about 10 5 to 10 6 
statistically independent configurations were generated. 
Fig. [1] shows typical snapshot pictures of such configura- 
tions for 4 choices of e. 



III. RESULTS 

When a single flexible polymer chain of chain length 
Nb tethered with one end to a flat impenetrable surface 
experiences an energy e that monomers win when they 
are adjacent to the surface, a transition of this "polymer 
mushroom" occurs from an essentially three-dimensional 
configuration to a quasi-two-dimensional "pancake" con- 
figuration when the adsorption energy exceeds a critical 
value e c [24 32] . This transition is a second order transi- 
tion, involving a nontrivial exponent <fr (the "crossover ex- 
ponent" |26[) such that the fraction of monomer-surface 
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FIG. 3. Ratio between the mean square gyration radius component perpendicular and parallel to the surface plotted 



the adsorption energy e for the bottle-brush backbone, R gbz /R gb << (a) and for the side chains, R 
monomer surface contacts plotted vs. adsorption energy e for the backbone of the bottle-brush, 
chains, N sc /(Nn c ); (d). 



vs. 

. 9C , z /i?g C ||(b). Fraction of 
N sb /N b (c) and for the side 



contacts N s /Nt, behaves as (in the limit Nt, — > oo) (2(| 

fcx \Nt\ e = e c (1) 

Despite a lot of effort (see 28 32] and references therein) 
the value of the (universal) exponent </> is not yet very 
precisely known, some estimates being compatible with 
4> = 0.50 ± 0.02 [H, & However, the most recent esti- 
mate from Monte Carlo simulations of the bond fluctu- 
ation model yielded [32[ 4> = 0.59, a value close to the 
early estimate (2(| tfi = 0.58. This adsorption transition 
also shows up in the chain linear dimensions, of course. 
While for e < e c all linear dimensions are of the same 
order [H, || 

R 2 gz oc R 2 gl{ oc N% v ,^« 0.588 , (2) 

for e > e c the perpendicular component of the gyration 
radius (related to the thickness of the "pancake" config- 
uration) remains finite [261 ] . 

Rl z <x(e/e c -l)-^«> , (3) 



while the parallel component exhibits a scaling compati- 
ble with two-dimensional self-avoiding walks [261 ] . 

R 2 gM oc (e/e c ~ i)^-")/* N ^ , V2 = 3/4 . (4) 

Hence when one plots the ratio R 2 , Z /R 2 || vs. e for several 
finite large values of Nb, one should see a family of curves 
which exhibit an intersection point for e = e c . 

Motivated by Eqs. H]), © and flU, we plot the fraction 
of monomer-surface contacts N s /N to t versus e (Fig.^a)), 
iVtot being the total number of effective monomers in 
the bottle- brush polymer, as well as the ratio R 2 Z /R 2 ^ 
(Fig. Efb)). Inspection of these "raw data" (Figs. |2T3")) 
already indicates that the adsorption transition occurs 
near e c (N) « 1 .00 ± 0.05 (see Fig. SI a close view of the 
ratio R 2 Z /R 2 ^ in the vicinity of e = 1); if e c (N) depends 
at all on side chain length N, the dependence is rather 
weak. However, it is remarkable that N s /N to t increases 
for e > e c (N) clearly much faster for the linear chains 
(N = 0) than for the bottle-brush polymers; N s /N tot 
seems to converge to a limiting function for large N that 
is independent of N. In contrast, the ratio R 2 Z /R 2 ^ 
varies more steeply the larger N, and there is no univer- 
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FIG. 4. Same figures as shown in Fig. O^b) and Fig. EJa) but for the data in the range 0.9 < e < 1.1. For N = 12 and N = 24 
the backbone chain length iVj, = 67 is too short so the expected scaling behavior (which requires an intersection of the curves 
at e c ) is not yet seen. 
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FIG. 5. Rescaled mean square end-to-end distance of the side chains, perpendicular to the surface, Rec,z/N 2v , vs. the normalized 
position k/(n c — - 1) in the coordinate system along the backbone (the n c side chains are labeled k = 0, 1, . . . , n — 1, the side 
chain with k — is next to the grafted backbone monomer, while the side chain with k = n c — 1 is next to the backbone free 
end). Five choices of adsorption energy e = 0, 1.0, 1.25, 1.5 and 2.0 are included, as indicated, as well as three choices of the 
backbone chain length, Nb = 67, 131 and 259, respectively. Case (a) refers to N — 6 and case (b) to N = 24. 



sal behavior; in particular, in the non- adsorbing regime 
Rg Z /R 2 g || increases monotonically with N. This behav- 
ior is plausible, since the longer the side chains are the 
more the excluded volume interaction with the planar 
substrate surface is felt. 

Of course, it would be desirable to extract more pre- 
cise estimates of e c (N) from a scaling analysis of the data, 
similar to the scaling analysis done in Refs. (2(| [29|, Hi], 
l32j . Unfortunately, as explained in the Appendix, no sig- 
nificant gain in accuracy in comparison with the simple 
intersection method of Fig. [^b) could be gotten. An ob- 
vious question with respect to the adsorption transition 
of bottle-brush polymers is the clarification whether the 
behavior of the side chains and of the backbone is fully 
analogous, or whether some characteristic differences oc- 
cur. Thus Fig. [3] disentangles in both the fraction of 
adsorbed monomers and in the ratio Rg Z /R 2 g || the con- 



tributions of the backbone from the contributions of the 
side chains. One sees remarkable differences between the 
behavior of the side chains and the backbone: while the 
ratio R gc z /R^ c ^ is essentially constant for e < e c and 
then decreases only rather gradually (almost linearly), 
the ratio R^z/^-gb \ \ shows most of its decrease before 
e c ~ 1.0 is reached. On the other hand, for large N a sig- 
nificant rise of N s b/Nb only starts at about e«1.2> e c . 
Some aspects of this behavior can be attributed to the 
fact that in the adsorption of bottle-brushes it easily hap- 
pens that the monomer density in the layer adjacent to 
the adsorbing wall gets rather high, thus there is not 
enough empty space in this layer to allow for the adsorp- 
tion of all the side-chain monomers. 

Since some of the side chains near the backbone chain 
end which is grafted to the substrate surface necessarily 
are rather close to the adsorbing surface, while this is not 
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FIG. 6. Same as Fig.[5j but for the component parallel to the surface, R 2 C u/N 2v , plotted vs. k/(n c — 1). 



(a) 



cm 



0.4 



0.3 



0.2 



0.1 



(b) 



£ = 0.0 






£=1.0 










£=1.25 










£=1.5 

£ = 2.0 








N b =259 
N b =131 
N b = 67 


□ 
+ 

A 


N = 6 


> 







0.5 

k/(n c -1) 




1 



cm 



£ = 0.0 




A ^N b =259 □ 
^a N b = 67 a 



0.1 




£ = 2.0 



0.5 

k/(n c -1) 



£ = 1.0 




FIG. 7. Same as Fig. [5] but for the rescaled mean square gyration radius of the side chains, perpendicular to the surface, 
RlcJN 2v , for N = 6 (a) and N = 24 (b). 



true for the side chains in the vicinity of its free ends, it 
also is interesting to resolve the side-chain properties as 
a function of their position along the backbone (Fig. [5]) . 
One can see that the linear dimensions of the side chains 
near to the grafted backbone chain end are strongly re- 
duced for non-adsorbed chains, while an analogous effect 
near the free backbone end is much smaller. No such 
asymmetry can occur for free bottle-brush polymers in 
bulk solution, of course, since there both backbone chain 
ends are strictly equivalent. It is remarkable, however, 
that this asymmetry effect vanishes almost completely 
when the bottle-brush polymer gets adsorbed: then also 
the free chain end of the backbone is close to the sur- 
face, and it does not matter whether it actually would 
be grafted or not. A similar asymmetry is found, how- 
ever, when one considers the mean square end-to-end dis- 
tance of the side chains parallel to the surface (Fig. [5]) : 
of course, then the magnitude of R 2 c ^/N 2v strongly in- 
creases with e, unlike R 2 C Z /N 2v , due to the tendency of 



the adsorbed polymers to form "pancake" conformations, 
cf. Eqs. ((H), and therefore these ratios systematically 
increase with increasing N. 

Similar observations can be made when one studies 
the mean square gyration radius components (Figs. [3 
|5J|. Thus, if one would normalize the data for R 2 c ^ with 

N 2V2 with ^ 2 = 3/4 instead of N 2l/ with v = 0.588, as 
done in Fig. [6j one finds plateau values for e = 2.0 in the 
center of the bottle-brush (0.2 < k/(n c - 1) < 0.8) at 
about 5.1 for N = 6 and at about 6.8 for N = 24. These 
values are of the same order as seen in Fig. [5] for e < 1.0, 
where the side chains still take the configuration of three- 
dimensional coils. Of course, side chain lengths N < 48, 
as studied here, are too short to reveal the asymptotic 
behavior R 2 c 1 1 oc TV 2 " 2 in the regime e > e c clearly. 

A rather unexpected feature, however, emerges when 
we examine the variation of the radii with e (Fig. (5]). The 
side chain linear dimensions are almost independent of 
Nf, for e < 1. The adsorbtion process for e > 1, however, 
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seems to be an TV-dependent two-step process: first a 
rather rapid decrease of R 2 gc z (and a rather rapid increase 
of R 2 gc ||) occurs, and the behavior for all three values of 
N is rather similar; then, near e w 1.25, the behavior 
changes, and the further variation depends distinctly on 

24 even is non-monotonic. Near 

We 



N and for R 2 gc z and N 
e w 1.7 the three sets of curves intersect each other 



speculate that for e < 1.25 the adsorption is controlled 
by the side chains near the grafted end of the backbone 
while for e > 1.25 it is controlled by all the side chains 
and the backbone, but clearly this point needs further 
study. Of course, Rg C ^ z can be considered as a measure 
of the thickness of the "pancake" . 

Alternative information on this thickness can be ob- 
tained when one simply studies the monomer density pro- 
file P(z), which again can also be resolved distinguish- 
ing contributions from the backbone only, from the side 
chains only, or even from the free chain ends of the side 



chains only (Figs. H0lfT4|) . 

While for N = 6 the monomer distribution of the 
bottle-brush as a whole for the non-adsorbed bottle- 
brush mushroom as well as the distribution of the back- 
bone and the side chain monomers are very similar to 
each other and similar to corresponding data for sim- 
ple non-adsorbed polymer mushrooms, the distribution 
gets a more complicated shape the longer the side chain 
length TV becomes. The strong shift of the maximum 
of all distributions with increasing N to larger z clearly 
can be attributed to the increasing backbone stiffness 
and the corresponding increase of the linear dimensions 
of the bottle-brush. When now the adsorption energy 
is switched on, the distributions show little change for 
e < 0.5 (not shown), while near the adsorption transition 
a profound change of the character of the distribution has 
occurred (Fig.[T2|): all distributions now exhibit a second 
maximum near the surface, for large enough TV, and also 
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FIG. 10. Monomer density profiles of the bottle-brush polymers P(z) as a function of the distance z from the surface, for 
backbone chain length Nt = 67 at the non-adsorbing case (e = 0.0) and four side chain lengths, N = 6, 12, 24 and 48. Case 
(a) refers to all monomers of the bottle-brush, case (b) includes monomers of the backbone only, case (c) includes only the side 
chains, while case (d) includes only the free ends of the side chains. 



the main maximum occurs now much closer to the ad- 
sorbing surface (but its position z max strongly increases 
with N, presumably because the bottle-brush needs to 
avoid "crowding" of monomers near the adsorbing sur- 
face). In the adsorbed regime (Figs. HUI H~4l one observes 
layering type oscillations near the wall, similar to the be- 
havior of off-lattice models for fluids close to hard walls. 
While the backbone monomers still have a pronounced 
maximum in P(z) away from the surface, Fig. I13f b) the 
corresponding results for the side-chain monomers do not 
show this maximum any more. This observation corrob- 
orates our interpretation of the non-monotonic variation 
of Rg C z in Fig. [SJ namely that adsorption occurs in two 
steps - first the adsorption of the monomers near the 
grafted end of the backbone takes place, while the ad- 
sorption transition of the rest of the bottle-brush is not 
yet completed at e = 1.25 for the side chain lengths N 
studied here. For e = 1.5, on the other hand, all profiles 
show only very little dependence on the side chain length 
(Fig. HI]): The lengths of still non-adsorbed "tails" and 
"loops" [27j then are rather short, and therefore the side 
chain length no longer is important. 

From the profiles, one can easily estimate the average 



thickness (z) defined as (z) — J zP(z)dz/ J P(z)dz, and 

o o 

study its dependence on e (Fig. [IS]). Again the two-step 
character of adsorption for long enough side-chains (N = 
24) is very clearly seen. 

As a final point of this study, we address the scattering 
function S(q) where the orientation of the wave vector q 
has been averaged over. Hence 

1 M AA 

= 7^EE ex P^- (*W*)]> (5) 

j=0 k=0 

where Af is the total number of monomers, from which 
the scattering is considered, and the average (...) in- 
cludes both a statistical average over the conformations 
of the bottle-brush and a spherical average over the di- 
rection of q. A distinctive advantage of the simulation 
is that the scattering from the total bottle-brush, the 
scattering only from the backbone, or only from the side 
chains are easily accessible. 

Fig. llGl correspondinglv compares data for S(q) consid- 
ering the total scattering from all the monomers of the 
bottle-brush, comparing data for 4 choices of N and 3 
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FIG. 11. Same as Fig. [TO] but for N b = 131 and e = 0.0 



choices of Nb, as indicated, and including data for both 
e = 0.0 and e = 2.0. While the data for e = 0.0 closely 
resemble the scattering function of isolated bottle-brush 
polymers in dilute bulk solution, as studied previously, 
revealing a law q~ x l v for small N and intermediate q, 
and a rod-like behavior q -1 for large q, for the strongly 
adsorbed case (e = 2.0) for large q the rod-like behav- 
ior is not seen for N = 24 and N = 48, rather one 
finds in the decade 0.1 < q < 1.0 a behavior propor- 
tional to q~ 1 l v ' 2 = <7~ 4 / 3 : the side chains in the strongly 
adsorbed case behave like two-dimensional self-avoiding 
walks: since there so many more monomers occur in the 
side chains than in the backbone, the rod-like characteris- 
tics of the latter is only revealed when one focuses on the 
scattering from the backbone only (Fig. [T7|) . Then one 
can clearly see two crossovers, Sb(q) ~ Nb(l — q 2 (R 2 g b )/3) 
at small q of order (R 2 b )^ 1 crosses over to q~ 4 / 3 in d = 2 
and to q -1 / - 588 i n d = 3 dimensions, while near q w 0.1 
a crossover to Sb{q) oc q^ 1 occurs, both for adsorbed and 
for non-adsorbed bottle-brushes. 



IV. SUMMARY 



In this paper a Monte Carlo simulation study of the ad- 
sorption of bottle-brush mushrooms" (i.e. bottle-brush 
polymers with a backbone chain end grafted to a flat 
structureless impenetrable surface) has been presented, 
using the bond-fluctuation model and assuming very 
good solvent conditions. The same range of backbone 
chain lengths (Nb < 259) and side chain lengths (N < 48) 
as used in a previous study of the same model in dilute 
solution in the bulk [19( has been used, since evidence has 
been presented that this range is fully appropriate to al- 
low a comparison with experiment [19( ■ Every backbone 
monomer carries one side chain. 

Both the backbone chain and the side chains are as- 
sumed to be fully flexible, and a short range attractive 
energy e (putting temperature ksT — 1 throughout) is 
assumed that attracts both monomers of the side chains 
and of the backbone to the surface in the same way. We 
show that near e c (N) = 1 the mushrooms cross over from 
a three-dimensional configuration (for e < 1) to the ad- 
sorbed, quasi-two-dimensional configuration for e > 1. 
These transitions occur roughly at the same value of 
e c (N) irrespective of N, but the latter has a pronounced 
effect deeper in the adsorbed region, where chain linear 
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dimensions reveal a two-step adsorption process (Figs. |H1 
[T5)) . We tentatively associate the first step of the ad- 
sorption (near e c ) to the side chains near the grafted end 
of the backbone, while the backbone as a whole for e 
not much larger than e c exhibits only a very small frac- 
tion N s ,b/Nb of adsorbed monomers (Fig. EJc)). Thus, 
most of the backbone monomers still occur in the tail or 
in rather long loops stretching away from the adsorbing 
surface. Only at somewhat larger values of e (namely 
for e » 1.2 in our case) the "trains" [27| of consecutively 
adsorbed backbone monomers get longer and more fre- 
quent, and the tail and the loops get correspondingly 
shorter. This behavior must cause an interesting inter- 
play with the behavior of the side chains: a side chain 
grafted to a backbone monomer that is part of an ad- 
sorbed "train" behaves basically like a small mushroom, 
and hence such a side chain gets easily adsorbed. How- 
ever, side chains grafted to a backbone monomer that 
belongs to the tail or to a large loop cannot yet get ad- 
sorbed easily, simply because typically the grafting site 
of this side chain is not close enough to the adsorbing 
surface. However, when e increases and the loops in the 
backbone would get smaller, the grafting site gets close to 
the adsorbing surface, and then the side chain easily can 
adsorb (and if it is long enough, the side chain adsorp- 
tion will have a kind of feedback effect on the backbone, 



dragging the remaining part of the non-adsorbed loop to- 
wards the surface as well). Clearly, this picture is qual- 
itative and somewhat speculative, but it suggests that 
through a cooperative interplay of side chains and back- 
bone a much more intricate behavior is possible than for 
the adsorption of linear polymers. One can even specu- 
late that due to the strong crowding effects upon adsorp- 
tion, occurring especially for long side chains, a complete 
adsorption of the bottle-brush can not occur from a con- 
figuration where both ends are already adsorbed and a 
train is still present in the interior of the polymer. Rather 
the adsorbed part towards the free end of the backbone 
would have do desorb again first and then the complete 
adsorption proceeds in a zipper-like fashion starting from 
the grafted chain end. Of course, development of an an- 
alytical model that could describe data such as shown in 
Figs . [21151 151 and [T51 would be highly desirable. 

We have also verified that strongly adsorbed bottle- 
brush polymers exhibit lateral linear dimensions that 
scale with the exponent of two-dimensional self- 
avoiding walks, and in the scattering function of the to- 
tal bottle-brush this behavior is seen as well (Fig. IT1)|) . 
However, the expected rod-like behavior on intermediate 
length scales is only seen when the scattering from the 
backbone is isolated (Fig. [17)) . Our data also imply rather 
different structures of adsorbed vs. non-adsorbed bottle- 
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FIG. 13. Same as Fig. HOI but for N b = 131 and e = 1.25 



brush polymers making it difficult to reason from ad- 
sorbed structures on the behavior of free bottle-brushes 
in solution! We emphasize this simple point, as a warn- 
ing to premature interpretations of corresponding exper- 
iments. 

Our study reveals many details (bimodal behavior 
of the density distribution P(z) of the monomers as a 
function of distance from the surface near adsorption 
fFig. rnil) . characteristic effects due to the backbone chain 
ends (Figs. [S] - |5}, etc., which all reflect in some way 
the interplay between the enthalpy won by adsorbing 
monomers, and various entropic effects. More work is 
clearly required to clarify the reasons for these detailed 
observations. Also, it would be very interesting to clarify 
the nature of the cross-sectional structure of the adsorbed 
bottle-brush on a coarse-grained scale (in the bulk, the 
bottle-brush can be viewed as a more or less flexible cylin- 
der: should we view the cross section of an adsorbed 
bottle brush as a sphere cap rather than a sphere?) A 
further interesting aspect is the question of the persis- 
tence length of adsorbed bottle-brushes [U |42| , and its 
relation to the persistence length of bottle-brushes in the 
bulk. Other interesting issues concern the effect of sol- 
vent quality [43j and possible irreversibility effects in the 
adsorption of the side arms Our study only con- 

siders adsorption under very good solvent conditions. It 



also is a very interesting problem (and relevant for ex- 
periment) to consider adsorption under poor solvent con- 
ditions where much more compact "pancake" structures 
should result [43[ rather than the configurations obtained 
here, that resemble two dimensional self-avoiding walks. 
Finally, we mention the possibility that for strong ad- 
sorption (e ^> e c ) the left-right distribution of the side 
chains with respect to the backbone is quenched, i.e., 
random fluctuations in this distribution have not enough 
time to relax. Such effects are predicted to have interest- 
ing effects on the conformation of such strongly adsorbed 
bottle-brush polymers (44|. We intend to study some of 
these issues in forthcoming work. 
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FIG. 14. Same as Fig. HOI but for N b = 131 and e = 1.50 




FIG. 15. Plot of (z) (note the logarithmic scale) vs. e, for various choices of Nb and N (labeled as bN b sN in the figure). Case 
(a) refers to all monomers of the bottle-brush, case (b) refers to backbone monomers only. 



Appendix A: Scaling Analysis of the adsorption 
transition bottle-brush polymers 

In the previous work where the adsorption transition 
of flexible chains end-grafted by a chain end on an im- 
penetrable surface was studied by Monte Carlo meth- 
ods [26|, [2^, [3l|, [lH , the location of the transition in the 



limit of infinite number of bonds between the monomers, 
Nf, — > 00, typically has been extracted from a scaling 
analysis [26| . The statement of scaling is that quantities 
such as the fraction N s /Nb of adsorbed monomers do not 
depend on the two variables Nb, n = (e — e c )Ac in the 
most general way, but are essentially a function of a sim- 
ple scaling variable, C — K N^, apart from a power law 
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(a) (b) 




FIG. 16. Log-log plot of the scattering function S(q) of the whole bottle-brush polymers versus wave number q, for side chain 
lengths TV = 6 (a), TV = 12 (b), TV = 24 (c) and TV = 48 (d). Three choices of backbone chain lengths TV, = 67, 131, and 259 
are included. Data for e = 0.0 are denoted as fcTV^eOO while data for e — 2.0 are denoted as &TV,e20. Straight lines indicate 
the power laws for rods three-dimensional coils (q^ 1 ^ with v = 0.588) and two-dimensional coils (g~ 1//y2 ,V2 = 0.75), 

respectively. 



prefactor. Thus 

N s = N$F S (() , k y , TV b y oo . (Al) 

The power laws quoted in the main text {Eq. (JT|)} simply 
result from Eq. (IA1|) for ( —> — oo, ( = 0, and £ — > +oo. 
Similar scaling results hold for the linear dimensions of 
the chains. E. g. the mean square gyration radius compo- 
nents of the backbone of the bottle-brush perpendicular 
and parallel to the surface are 

R 2 gb , z = N^F bz (0 , R 2 gbM = N^F b]] (0 , (A2) 

and again Eq. (|A2[) is supposed to be valid in the limit 
where both k — > and TV b — > oo. For Eq. (|A2|) . the 
scaling functions for £ < are less interesting, since 
F bz (C -> -co), F b ||(C -> -co), F 6 ,(0), and F b]] (0) are 
finite constants {cf. Eq. (|2])}. However, for £ — > +oo one 
has 

^(C) oc C 2 ^ , Fn| (0 « C 2 ^ 2 -^/* , (A3) 

and combining Eqs. (IA2[) . (|A3I) one hence recovers 
Eqs. ©, (0]) of the main text. 



Of course, the singularities described by Eqs. <HJ)- 
of the main text and the above Eqs. (IA1|) - ()A3|) emerge 
only in the double limit n — > 0, TV;, — > oo: for any fi- 
nite TV;, functions TV s (e), R 2 gb 2 (e) and R 2 gb ,,(e) as regular 
functions of e, the singular behavior is rounded off, as is 
well known 26] . Thus the estimation of e c is a nontrivial 
matter. What usually is done, see e.g. Refs. (26ll3l1, |32]| . 
is a "data collapsing method" : one varies both e c and the 
(originally unknown) value of the crossover exponent to 
obtain an optimal "data collapse" of the set of functions 
TV s (e,TVb), R 2 gbtZ {e,N b ) and R 2 b ^ (e, N b ) on these "master 
curves" which represent then the three scaling functions 
F S (C), F bz (C), and F 6N (C). Figs. HMD show some typical 
attempts to do this with our data (each figure has four 
parts to show the results for the four side chain lengths 
TV = (no side chains), TV = 6, TV = 12, and TV = 24, re- 
spectively). Here we have used as a further constraint the 
"universality principle" , namely the crossover exponent 
<t> should be a universal constant, independent of irrele- 
vant microscopic details such as the length TV of the side 
chains (as long as TV is finite, while the limit N b — ► co 
is considered). Note also that the ratio R 2 gbz /R 2 b i ac- 
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(a) (b) 




FIG. 17. Log- log plot of the scattering function Sb (q) of the monomers in the backbone vs. wave number q, for side chain 
lengths N = 6 (a), N = 12 (b), TV = 24 (c) and TV = 48 (d). Three choices of backbone chain length N b = 67, 131, and 259 are 
included in the same notation as in Fig. [16] is used. 



cording to Eq. (|A2[) should be a function of £ alone, and 
plotted vs. e should have a unique intersection point at 
e c . (However, as shown in Fig. H] one would need data 
for larger Nh and very good statistical accuracy to show 
this). 

The straight lines included in these log-log plots of the 
resulting estimates for the scaling functions that they 
must converge to for | £ |— > oo yield the power laws 
quoted in Eqs. Since v = 0.588 and v 2 = 3/4 

is independently known, it is again only the exponent <f> 
(as in the scaling variable kN?) which matters. To sim- 
plify matters, we only present the two extreme choices 
for <j) here: the estimate <j> = 0.5 (which happens to coin- 
cide with the crossover exponent for adsorption of Gaus- 
sian chains [26I428I . but is supported by a renormalization 
group estimation [30| and some Monte Carlo studies [291 ]) 
and the most recent estimate <f> = 0.59 from Monte Carlo 
studies [32[ . Comparison of Figs [TS1 US] for N s shows 
that we confirm the finding of Descas et. al. 32], that 
(f> = 0.59 yields a slightly better data collapse on the mas- 
ter curve F s ((), e c (N) being always chosen such that for 
e < e c the master curve has the correct slope. For e > e c 
the data for large C fall systematically somewhat below 
the expected power law, if <f> = 0.5 is chosen. However, 



as a caveat we mention that for finite Nh the data are 
not expected to follow the master curve for very large £, 
since N s /Nh — > 1 (for the case N — and large enough e, 
where the chain is adsorbed in a two dimensional config- 
uration). So, it must happen that for large £ the curves 
bend away from the power law that describe the asymp- 
totic limit of the scaling function, and this "saturation 
effect" occurs the later the larger Nh is, and this is ex- 
actly what one sees in Fig. [hSr a). while in Fig. [19] (a) the 
same effect is somewhat less pronounced. However, in 
our opinion the evidence from these plots that tj> = 0.59 
is "better" than <f> = 0.50 is somewhat weak, and much 
longer chains would be needed to reach a really firm con- 
clusions. 

When one chooses <fr = 0.50, the resulting estimates 
for e c (N) would be 0.98, 1.01, 1.02, and 1.00, for side 
chain lengths N = 0, 6, 12, 24, respectively. When one 
chooses 4> — 0.59, on the other hand, most of the esti- 
mate are a little bit larger, namely 1.00, 1.03, 1.05, and 
1.00, respectively. As expected, the best estimates for <f) 
and e c (2V) are correlated. However, it must be admitted 
that slightly different choices of e c (N) than those that 
are shown here yield a data collapse that is only slightly 
worse. Given the fact that Eqs. ([I])-® and (|A"l J) -(|A"3" l) 



15 




FIG. 18. Log-log plot of the rescaled surface contacts N s /N{ ot vs. | k \ N{ ot for cj> = 0.50. Data for backbone lengths Nb = 67, 
131, 259, and 515 are included, as well as side chain lengths (a) N = 0, (b) N = 6, (c) N = 12, and (d) N = 24. 



are only asymptotically valid in the double limit K — > 
and Nf, — > oo, while for finite Nb and finite nonzero k cor- 
rections to scaling may be present, a reliable judgment 
of accuracy for the estimates e c (iV) is difficult. For this 
reason, we have quoted e c « 1.0 as an estimate for all 
choices of N studied in the main text. 

In fact, the estimation is not really improved when one 
considers the scaling of the mean square gyration radius 
(Figs. [201 ED- For N = and N = 6, the quality of the 
data collapse for <fi — 0.50 and 4> = 0-59 is of comparable 
quality. For N = 12 and N = 24, scaling seems to work 



well only for large £, not for small C tl indicating clearly a 
more complex behavior right at the adsorption transition. 
The data show that the choices for e c (N) quoted above 
are compatible with the behavior of the chain dimensions 
too, at least roughly. 

The result that the addition of side chains (that are 
still short, of course) has so little effect on the location 
of the adsorption transition clearly is surprising and un- 
expected. Therefore we have presented in this Appendix 
our evidence for this fact in some detail, not withstand- 
ing the uncertainty about the best value of the crossover 
exponent; this problem is not solved by our analysis. 



[1] M. Zhang and A. H. E. Miiller, J. Polym. Sci. Part A: 

Polym. Chem. 43, 3461 (2005). 
[2] A. V. Subbotin and A. N. Semenov, Polymer Science, 

Ser. A 49, 1328 (2007). 
[3] S. S. Sheiko, B. S. Sumerlin, and K. Matyjaszewski, 

Progr. Polym. Sci. 33, 759 (2008). 
[4] I. I. Potemkin and V. V. Palyulin, Polymer Science, Ser. 

A 51, 123 (2009). 
[5] Y. Tsukahara, K. Mizuno, A. Segawa, and Y. Yamashita, 

Macromolecules 22, 1546 (1989); Y. Tsukahara, K. Tsut- 



sumi, Y. Yamashita, and S. Shimada, Macromolecules 

23, 5201 (1990). 
[6] M. Wintermantel, M. Gerle, K. Fischer, M. Schmidt, I. 

Wataoka, H. Urakawa, K. Kajiwara, and Y. Tsukahara, 

Macromolecules 29, 978 (1996). 
[7] R. Djalali, S.-Y. Li, and M. Schmidt, Macromolecules 35, 

4282 (2002). 

[8] T. Stephan, S. Muth and M. Schmidt, Macromolecules 

35, 9857 (2002). 
[9] C. Li, N. Gunari, K. Fischer, A. Janshoff, and M. 



16 




FIG. 19. Log-log plot of the rescaled surface contacts N s /N{ ot vs. | k \ N{ ot for cj> = 0.59. Data for backbone lengths Nb = 67, 
131, 259, and 515 are included, as well as side chain lengths (a) N = 0, (b) N = 6, (c) N = 12, and (d) N = 24. 



Schmidt, Angew. Chem. Int. Ed. 43, 1101 (2004). 
[10] R. V. Iozzo (ed.) Proteoglycans: Structure, Biology, and 

Molecular Interactions (Marcel Dekker, New York, 2000). 
[11] N. C. Kaneider, S. Dunzendorfer, and C. J. Wiedermann, 

Biochemistry 43, 237 (2004). 
[12] J. Klein, Science 323, 47 (2009). 

[13] A. Sartori, A. Johner, J.-L. Viovy, and J.-F. Joanny, 

Macromolecules 38, 3432 (2005). 
[14] P. G. Khalatur, A. R. Khokhlov, S. A. Prokhorova, S. S. 

Sheiko, M. Moller, P. Reineker, D. G. Shirvanyanz, and 

N, Starovoitova, Eur. Phys. J. E. 1, 99 (2000). 
[15] I. I. Potemkin, A. R. Khokhlov, S. Prokhorova, S. S. 

Sheiko, M. Moller, K. L. Beers, and K. Matyjaszewski, 

Macromolecules 37, 3918 (2004). 
[16] N. Gunari, M. Schmidt, and A. Janshoff, Macromolecules 

39, 2219 (2006). 
[17] S. Rathgeber, T. Pakula, A. WilK, K. Matyjaszewski, 

and K. L. Beers, J. Chem. Phys. 122, 124904 (2005). 
[18] B. Zhang, F. Crohn, J. S. Pedersen, K. Fischer and M. 

Schmidt, Macromolecules 39, 8440 (2006). 
[19] H.-P. Hsu, W. Paul, S. Rathgeber, and K. Binder, Macro- 
molecules 43, 1592 (2010). 
[20] I. I. Potemkin, Macromolecules 39, 7178 (2006). 
[21] S. Elli, G. Raffaini, F. Ganazzoli, E. G. Timoshenko and 

Y. A. Kuznetsov, Polymer 49, 1716 (2008). 
[22] P. G. de Gennes, Scaling Concepts in Polymer Physics 

(Cornell University Press, Ithaca, 1979). 



[23] Note that previous work on the adsorption transition of 

comb polymers used N = 5 only [2l| . 
[24] P. G. de Gennes, J. Phys. (Paris) 37, 1445 (1976). 
[25] P. G. de Gennes, Macromolecules 13, 1069 (1980). 
[26] E. Eisenriegler, K. Kremer, and K. Binder, J. Chem. 

Phys. 77, 6296 (1982). 
[27] G. J. Fleer, M. A. Cohen-Stuart, J. M. H. M. Scheut- 

jens, T. Cosgrove and B. Vincent, Polymers at Interfaces 

(Chapman and Hall, London, 1993). 
[28] E. Eisenriegler, Polymers Near Surfaces World Scientific, 

Singapore, 1993)/ 
[29] R. Hegger and P. Grassberger, J. Phys. A: Math. Gen. 

27, 4069 (1994). 
[30] H. W. Diehl and M. Shpot, Nucl. Phys. B528, 595 (1998). 
[31] S. Metzger, M. Miiller, K. Binder and J. Baschnagel, 

Macromol. Theory Simul. 11, 985 (2002); J. Chem. Phys. 

118, 8489 (2003). 
[32] R. Descas, J.-U. Sommer, A. Blumen, Macromol. Theory 

Simul. 17, 429 (2008); J. Chem. Phys. 120, 8831 (2004). 
[33] V. A. Ivanov, J. A. Martemyanova, M. Miiller, W. Paul 

and K. Binder, J. Phys. Chem. B113, 3653 (2009). 
[34] I. Carmesin and K. Kremer, Macromolecules 21, 2819 

(1988). 

[35] H. P. Deutsch and K. Binder, J. Chem. Phys. 94, 2294 
(1991). 

[36] W. Paul, K. Binder, D. W. Heermann and K. Kremer, J. 
Phys. II (Paris) 1, 37 (1991). 



17 



(a) 



10 



&x> 1 



10" 



s - 1 0" 2 
cc 



10"' 







II 

- + □ X 


slope = 2(v 2 -v)/<|) 


+ □ * A 

z 


* 


b67s0 


+ [i slope = -2v/(() 


b131s0 




. b259s0 




b515s0 


e c =0.98, i)>=0.5o\^ 



10" 



10 



kN2 



10" 



(b) 



10" 
10 

1 



z 10 



10 z 



10"' 
10"' 



II 

+ □ *+ 


iXribSflW** 8 *"*" 

^ * A sk>pe = 2(v 2 -v)/<|> 


z - □ x + 


*b + 

^^itoSlope = -2v/<|) . 




b67s6 + 




" b131s6 




_ b259s6 * 


e c =1.01, <h0.50^N. 



10" 



10 



kN2 



(c) 



10" 
10 

z 1 

? 10" 1 

Z -2 

"5 10^ 

DC -q 

10 d 



10 



10 





slope = 2(v 2 -v)/(|) 


+ 




z 


*° ^fx^lope = -2vA|> 


b67s12 




- b131s12 


x e c =1 . 02 ■^a^S'^ 


. b259s12 




)" 1 


1 10 







(d) 



10" 
10 

z 1 

? 10" 1 

Z -2 

"5 10^ 
CC -q 

10 d 



10 



slope = 2(v 2 -v)/(|) 



2 x □ □ + n Q 

* * * * slope = -2v/«|) 

+ x 



b67s24 + 
b131s24 □ 

b259s24 x e c =1 .00; <|)=0.50 




10 



1 



10 



FIG. 20. Log-log plot of the mean square gyration radius components perpendicular and parallel to the surface, B%b !Z /Nb" and 
B? gb ||/jVj", respectively, vs. nNf ot for <j> — 0.50. Data for backbone lengths Nt = 67, 131, 259, and 515 are included, as well as 
side chain lengths (a) TV = 0, (b) TV = 6, (c) TV = 12, and (d) TV = 24. Here the Flory exponents v 2 = 3/4 (2D), and v = 0.588 
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FIG. 21. Log-log plot of the mean square gyration radius components perpendicular and parallel to the surface, B%b !Z /Nb" and 
B? gb ||/jVj", respectively, vs. nNf ot for <j> — 0.59. Data for backbone lengths TV, = 67, 131, 259, and 515 are included, as well as 
side chain lengths (a) TV = 0, (b) TV = 6, (c) TV = 12, and (d) TV = 24. Here the Flory exponents v 2 = 3/4 (2D), and v = 0.588 
(3D). 



